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(~) Abstract 

Wc examine the effects of the different viscosity prescriptions and the magnitude of the viscosity pa- 
rameter, a, on the structure of the slim disk, and discuss the observational implications on accretion-flow 
structure into a stellar-mass black hole. In contrast with a standard disk, in which the "a" value does 
not affect significantly the local flux, radiation from the slim disk is influenced by the a value. For the 
range of a = 10~ 2 ~ 10° we calculate the disk spectra and from spectral fitting we derive TJ° bs , maxi- 
mum temperature of the disk, i?f bs , the size of the region emitting blackbody radiation with T]° b8 , and 
p = — dlnT ff/dln r, the slope of the effective temperature distribution. It was founded that the estimated 
7j° b8 slightly increases as a increases. This is because the larger the magnitude of viscosity is, the larger 
' becomes the accretion velocity and, hence, the more enhanced becomes advective energy transport, which 

means less efficient radiative cooling and thus higher temperatures. 

Furthermore we check different viscosity prescriptions with the form of the viscous stress tensor of 
Ptotai, where (3 is the ratio of gas pressure (p gas ) to total pressure, ptotai (=P g as +Prad), and fi is 
a parameter (0 < [i < 1). For \i = we have previously found that as luminosity approaches the Eddington 
, Le,R°^ s decreases below 3r g (3r g corresponds to the radius of the marginally stable circular orbit, r ms , 
with r g being Schwarzschild radius) and the effective temperature profile becomes flatter, T c ff oc r -1 / 2 . 
Such a slim-disk nature does not appear when \x is large, \i ~ 0.5, even at the Eddington luminosity. Hence, 
the temperature of the innermost region of the disk sensitively depends on the /i value. We can rule out the 
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case with large fx (~ 0.5), since it will not be able to produce a drop in Rf^ with an increase in luminosity 
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as was observed in an ultraluminous X-ray source, IC 342, source 1. 
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1. Introduction discussed in relation to their X-ray properties since late 

1980's (Fabbiano, 1989). For instance, Ultra-Luminous 

The accretion disk model by Shakura and Sunyaev Compact X-ray Sources (hereafter, ULXs) in nearby spi- 

(1973, hereafter, SSD) has been regarded as the "stan- ral galaxies are very mysterious objects. Their luminosi- 

dard" model for the low mass X-ray Binaries (LMXBs) ties amount to around 10 39_40 crg s _1 , which means the 

and galactic black-hole candidates (GBHCs) in their soft mass (M) of black holes, if they exist, should exceed 10- 

state. Their analytical solutions are self-consistent under 100 M Q . However, the MCD (multi-color disk, Mitsuda 

some specific assumptions. However, when its luminos- et al. 1984) fitting derives relatively high temperatures, 

ity L approaches the Eddington luminosity, Le, some as- T°^ s ~ 1.0 — 2.0 keV (Mizuno et al. 1999, Makishima 

sumptions break down and thus the standard model can- et al. 2000), compared with the theoretical expectations 

not adequately describe the nature of such high luminosity based on the SSD model, which is significantly below 1.0 

objects. The critical luminosity is L ~ 0.3Le (Abramowicz keV for M > 100M Q . Here, T]° b8 is the color tempcra- 

ct al. 1988, Laor, Netzer 1989). This fact is theoretically ture derived from the observed spectrum (for the exact 

established but not widely known and still most of in- definition, sec Section 3.3.). In our Galaxy, there ex- 

vestigations have been conducted by simply applying the ists a famous galactic microquasar, GRS 1915+105, which 

canonical SSD even to near-Eddington luminosity systems shows the complicated time/spectral variation. Its lumi- 

without any comments. nosity is near the Eddington luminosity in its high state 

Recent observations have detected several kinds of ex- > a few x 10 39 erg s _1 (Belloni et al. 2000). These ob- 

otic supcr-Eddington sources, whose luminosity seemingly jects can never be explained in the framework of the SSD 

exceeds the Eddington luminosity. The properties of their theory. 

energy source cannot be described by the SSD theory, On the theoretical side, there has been growing recog- 

and hence its origin is still an open astrophysical qucs- nition that the super-critical accretion disk model, so- 

tion. The nature of super-Eddington sources has been called "slim disk" (Abramowicz et al. 1988) may be a 
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correct model to describe these high luminosity objects. 
Its significance has been re-recognized in connection with 
luminous AGNs (Szuszkicwicz et al. 1996, Wang et al. 
1999), narrow-line Seyfertl galaxies (NLSls; Mineshigc ct 
al. 2000), and GBHCs shining near the Eddington lumi- 
nosity (Fukue 2000, Watarai et al. 2000), since the disk 
is stable due to the advective entropy transport, even at 
the super-critical accretion rates. The ULXs are discussed 
along this line in the first time by Watarai et al. (2001) 
and it has been shown to explain the basic observed ten- 
dency, although more detailed analysis is needed (Watarai 
in preparation). 

The local stability against thermal perturbations can 
be discussed in terms of the thermal equilibrium curve, so 
called "S-shaped" curve (Taam, Lin 1986, Abramowicz et 
al. 1995). The local stability properties for the station- 
ary slim disk models with different viscosity prescriptions 
were discussed by Szuszkicwicz (1990), who found that 
when "p" increases the thermally unstable region in the 
disk becomes smaller. Here, we prescribe viscous-shear 
tensor, the different form of shear stress tensor to be, 
t rtp = -ap^asPtotai' is the S as pressure, and p to tai is 

the total pressure), where p is the parameter (0 < p < 1). 

The global stability analysis for slim disk model has 
been carried out by Honma et al. (1991) who found that 
the disk was stabilized as the viscosity parameter "/i" in- 
creased. Szuszkiewicz and Miller (1997) confirmed this 
finding, showing that the models are stable not only when 
they are locally stable but also when the local analysis 
predicts an unstable region with radial dimension smaller 
than the shortest wavelength of the unstable modes. In 
the most recent study by Zampieri et al. (2001) have cal- 
culated the spectra of a slim disk model along one full 
thermal limit cycle. However, there is little research that 
carried out direct comparison with X-ray observations. 

Our concern in the present study is how the viscosity 
prescription and its strength act on the structure of ac- 
cretion disks, and how they manifest themselves in the 
fitting parameters (7]° bs , i£° bs ) derived from the observa- 
tional spectra in soft X-ray band. If T|° bs and i?? bs can be 
accurately determined from the observation, conversely, 
the functional form of disk viscosity might be able to con- 
strain. A comprehensive understanding of the structure 
and spectrum of slim disks is an outstanding issue. 

In the next section we introduce the basic equations and 
explain our calculation method. In section 3, the results 
of our calculation are presented for a variety of viscos- 
ity parameters (a,p). The observational features of the 
slim disk arc summarized by using the fitting parameters, 
T|° bs , i?°, bs , and p in section 4. Final section is devoted to 
discussion and concluding remarks. 

2. Basic Equations 

Basic equations we use here are all the same as those 
in the previous paper (Watarai et al. 2000) except for the 
viscosity prescription (see Kato et al. 1998, Chapter 8 for 
the detailed description). We solved the height-integrated 
equations (Hoshi, 1977) in the radial direction on the basis 
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of pseudo-Newtonian potential (Paczyhsky, Wiita, 1980), 
and neglected self-gravity. The pseudo-Newtonian poten- 
tial is not an excellent approximation around the event 
horizon ~ r g . Thus, caution is needed when the derived 
is small, R°^ s < r g . Other general relativistic effects 
needs also to take account, then. Otherwise, our approx- 
imation does not introduce serious errors. 

Honma et al. (1991) were the first to carry out the 
calculation above circumstances (including the effect of 
viscosity prescriptions of the slim disk). The main differ- 
ence of our calculation and theirs is that they truncated 
the disk inner edge at 2.7 r g , while we integrate the disk 
down to the real vicinity of the event horizon at r g . 

The momentum equation in the radial direction is 

dv r 1 dW l 2 -l\ W din fl K 

v ^ + ^ = — Y~d^~ 7 (1) 

where v r is the radial velocity of gas, £ = 2Ij,p^H is surface 
density, and W = J ptotaidz = 2I^pqH is pressure integral, 
H is disk scale height, respectively. Here, the pressure, 
P0j and the density, p$, are related to each other by the 
relation p oc Pq +1 ^ N (the subscript '0' indicates the quan- 
tities on the equatorial plane). We fixed the polytropic 
index to be N = 3, for which the numerical constants are 
7 3 =16/35 and 7 4 =128/315 (Hoshi, 1977). 

The specific angular momentum and the Keplerian an- 
gular momentum are defined by £ (= rv^) and £k = ?~ 2 ^k, 
respectively, where fix is the Keplerian angular frequency. 
The angular momentum conservation is 

M(£-£ ia ) = -2nr 2 T rv , (2) 

where M is mass accretion rate, 

M = —l-KTVrYi = constant., (3) 

from continuity equation and T rip is the height integrated 
viscous stress tensor, T rip = J t rv dz. 

Hydrostatic balance in the vertical direction leads 

H 2 nl = (2N + 3)^ = cl (4) 

where c s is the sound speed. 
The equation of state is 

Ptotai = Pgas + Prad = — P0 T + ^ Tq , ( 5 ) 

P O 

where R, p, and a are the gas constant, the mean molec- 
ular weight, and the radiation constant respectively. As 
to the viscosity prescription we rewrite it in terms of the 
effective a; 

T rip = -a e sW with a off = = a(p gas /ptotai) M - (6) 

This form is convenient to compare with the cases with 
the usual viscosity prescription (i.e., p = 0). If we take 
an extreme case of p — » limit, the stress tensor is pro- 
portional to total pressure; i.e., the same as that of the 
SSD prescription. In the large p limit (p — »• 1), on the 
other hand, the shear stress tensor depends solely on gas 
pressure. This situation may realize when there grows tur- 
bulence generated by chaotic magnetic fields (Lightman, 
Eardley 1974; Sakimoto, Coroniti 1981). 
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The energy equation involves the viscous heating, ra- 
diative cooling and advective cooling terms; namely, 
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Here, sq is the entropy on the equatorial plane, k is 
the Rosseland-mcan opacity, R = n cs + Kg = 0.40 + 0.64 x 
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g cm with p and T being averaged den- 



sity and temperature, respectively, (p= 16/35po and T 
2/3Xb) and other symbols have their usual meanings. We 
solved the above set of equations (||)-(^|). The solution is 
calculated from the outer edge of the disk (r out = 10 4 r g ) to 
the vicinity of the central black hole through the transonic 
point. We adjust l m , specific angular momentum which 
is finally swallowed by the black hole, so as to satisfy the 
regularity condition at the transonic radius. 

For evaluating the effective temperature, we assume 
that all the radiative energy is emitted with blackbody 
as a first approximation. Namely, using the Stefan- 
Boltzmann's law, Q~ ad = 2F = 2oT 4 ff , we define the lo- 
cal effective temperature to be T e g = (Qrad/2cr) ' 25 . As 
is widely known, however, the effect of Compton scatter- 
ing is substantial in electron scattering dominated layers 
(Czerny, Elvis 1987; Ross et al. 1992; Shimura, Takahara 
1995; Wang et al. 1999). Wc include this effect by in- 
troducing a spectral hardening factor (see Section. 3.3.). 
Further, photon trapping effect (e.g., Mineshige et al. 
2000) needs to be considered. However, we ignore these 
effects in order to understanding the basic tendency of the 
spectrum. 

Finally, wc define the normalized parameter m = 
M/Mcrit, with M crit = L E /c 2 = 1.3 x 10 17 (A//M Q )g s" 1 . 
This definition may differ in different literatures. In our 
notation, disk luminosity corresponds to the Eddington 
luminosity, when the accretion rate is m = 16 (provided 
that the pseudo-Newtonian potential is used). We vary 
the normalized parameters, m = M/Mq, m, a, and p, and 
calculate the structure of optically thick accretion disk for 
a wide range of parameters. 

3. The effects of varying viscosity prescriptions 

3.1. The a-Dependence 

We first examine the effects of changing the a parame- 
ter within the SSD prescription, t rip = — aptotai, i-e., p = 0. 
It is well known that a does not affect significantly the 
spectrum in the standard disk theory, since total radia- 
tion output rate is balanced with the energy generation 
rate via viscous processes and does not explicitly depend 
on density nor velocity. In other words, effective temper- 
atures can be nearly uniquely determined by black-hole 
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Fig. 1. The effective temperature distributions and the 
emergent spectra for different a parameters (cv=0.01, 0.1, and 
1.0). Solid lines represent the calculated models (for m = 10 
and 100). Dashed lines are the results based on the stan- 
dard-disk relation (Shakura, Sunyaev 1973). For a small ac- 
cretion rate, m = 10, the spectrum does not largely depend 
on a parameter, where as at a high accretion rate, m = 100, 
a dependence appears in the temperature profile in the inner 
region. Likewise, the spectral changes are not appreciable in 
the case of m = 10 while in m = 100 case, the spectrum has 
a dependence reflecting the a dependence of the temperature 
distribution. Black- hole mass is fixed to be IOMq. 

mass and mass-accretion rate as functions of radius. In a 
high accretion-rate system, in contrast, a dependence ap- 
pears, since radiative cooling is no longer balanced solely 
with viscous heating so that the disk inner temperature 
can change even for the same M and M. 

Figure |l| shows the temperature profiles and spectra for 
different a and m. In the temperature profiles, we found 
T oc r~ 3 / 4 at low m ^ 10 (called as the SSD regime) and 
T oc r~ x l 2 at large to > 100 (called as the slim-disk regime, 
see Watarai, Fukue, 1999, Watarai et al. 2000; Wang et 
al. 1999). For instance, when a increases by one order 
of magnitude, the temperature in the vicinity of the disk 
inner region only slightly rises by a factor of ~ 2, and 
X-ray radiation spectra get a bit harder. 

Next, we check the relative importance of advective en- 
ergy transport in the energy equation. Figure || plots the 
ratio, Qadv/Qvis' f° r different a and to. When a increases, 
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Fig. 2. The ratio of advectivc energy transport to the vis- 
cous heating rate for different accretion rate (m=33 and 
100). For all the cases, the viscous heating rate Q^ la is bal- 
anced with the radiative cooling Q~ ad at larger radii, thus 
|Q — . | -C Q^ is . At higher accretion rates, the advective cool- 
ing Q~ dv is substantial at inner radii, but still, the ratio is 
mostly less than 0.5;that is, the flow is not advection domi- 
nated. Black-hole mass is fixed to be M = IOMq , and we set 
a = 0.1 (the solid lines) and 0.01 (the dotted lines). 

this ratio also increases; that is, the large a is, the more 
effectively advection cooling works. We note, however, 
that even for rh = 100, still disk is not entirely advection 
dominated (the maximum ratio is ~ 0.5). Thus, it is not 
appropriate to call the slim disk the optically thick ADAF 
(see also Abramowicz et al. 1988). 

Figure 1 explicitly demonstrates that the disk inner 
edge decreases as rh increases. Why, then, can the in- 
ner edge of the disk be smaller than the radius of the 
marginally stable circular orbit, r ms ~ 3r g ? One of the 
reasons is a decrease of the transonic radius of the flow, rs , 
( r s r ms) with increase of accretion rate; i.e. rs / r g =3.00, 
2.84 and 2.69 for m < 10,33 and 100, respectively. We ob- 
tained rs for several rh and confirmed that the tendency of 
rs is consistent with that of Abramowicz et al. (1988). In 
the case of super-critical accretion, moreover, a pressure 
gradient force is important as in the case of a thick torus 
(Abramowicz et al. 1978). For small a (< 0.1) accretion 
is, in fact, driven by pressure gradient. Then, the pressure 
gradient force tend to push material near the inner edge 
of the disk inward to smaller radii, even if the material 
has larger angular momentum than the equilibrium value, 
for which centrifugal force is balanced with gravity force. 
Finally, even when the flow is unstable for circular mo- 
tion, the growth time is on the order of free-fall accretion 
time, while in the meantime mass is continuously supplied 
from outside with large (comparable to sound) velocity at 
large rh. There always exists substantial material inside 
r ms , emitting radiation. 

For large a value (~ 1.0) the accretion is driven by the 
viscosity, however, such a flow can be effectively optically 
thin. It was pointed out by Beloborodov (1998) that in 
the case of very large a ~ 1.0, the inner region of the disk 
may be overheated due to decreasing the effective optical 
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Fig. 3. The same as figure 1 but different fi parameters (fi 
= 0, 0.1, 0.4, and 0.49). 

depth, but strong convective motion will reduce thermal- 
ization timcscale, which may avoid overheated problem 
for moderately large a ~ 0.1 (Mineshige et al. 2000). 

3.2. The fj. Dependence 

In this subsection, we discuss the fi dependence of the 
slim disk structure. In equation (^), we fix a (=0.1) and 
M (= 10M Q ) but change the value of /j, (=0, 0.1, 0.4, and 
0.49). If we take the limit of fi approaching unity; that is 
t rip = — apg as , and then the disk will be completely stabi- 
lized, since then temperature dependence of the heating 
rate (oc t rtp ) is now less than that of the cooling term, 
a condition for thermal stability (see Kato et al. 1998, 
Chap. 4). 

Figure ^ display the temperature profiles and the spec- 
tra for different \i values. When m = 10, the temperature 
does not show \i dependence. In the rh = 100 case, in 
contrast, the temperature appreciably decreases around 
transonic region (< 3r g ) as fi decreases, keeping the same 
temperature slope, T oc r -1 / 2 , at large radii, r ^> 3r g . 

The basic interpretations of these results are that a large 
ii tends to decrease a e g at smaller radii, thus decelerating 
infall motion of gas into black hole. This rather promotes 
radiative cooling, thereby decreasing temperature. This 
means, the parameter \i more appreciably affects the flow 
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structure in the vicinity of the central hole. 

3.3. Observable fi effects 

Finally, we show the effects of changing \x more explic- 
itly in the model spectral fitting diagram in comparison 
with the observational properties of soft X-ray dominated 
sources. The method of model fitting is the same as that 
described in Mineshige et al. (1994; see also Watarai et al. 
2000). The fitting parameters are i?|^ od , which determines 
the height of the spectral peak, X|™ od , specifying the fre- 
quency of the spectral peak, and the power-law index in 
the temperature distribution, "p"(= —dlnT/dlnr), which 
controls the spectral slope at low-energy (the Rayleigh- 
Jeans) side of the peak. We thus assume the form, 

where T™° d is the temperature at the outer edge and 
B(E,T) is the Planck function. The limit of p ~ 0.75 
corresponds to the original MCD model, while p ~ 0.5 in- 
dicates the typical temperature gradient of the slim disk 
model, T oc r -1 / 2 . The results of the fitting for different 
\i values are summarized in Tables 1-3 (see also Tables 
4 and 5 for different a cases). Throughout this section, 
we fix a = 0.1. Note that we are more concerned with 
relative changes of the fitting parameters as functions of 
M and M, rather than their absolute values for specific 
values of M and M, since realistically, general relativis- 
tic effects and other effects (such as self-shielding and the 
disk tilting) should modify the results, which are not easy 
to evaluate precisely. However, general tendency does not 
change as long as the central black hole is non-rotating 
(see Manmoto, Mineshige 2001 for the detailed discus- 
sion). 

According to our previous paper (Watarai et al. 2001), 
the "X|° bs — i?? bs " diagram is useful for diagnosing the ob- 
jects having a soft thermal component in the X-ray spec- 
trum; with this we can constrain various models (see also 
Belloni 2000, Mizuno et al. 2000). We considered the 
spectral hardening factor, £ (Shimura, Takahara 1995), 
and the correlation factor, r\ (Kubota at al. 1998), in the 
fitting results. The hardening factor £ represents the spec- 
tral hardening due to Comptonization processes within 
the disk (i;° bs = £T™ od , see the second column of tables). 
The rj parameter is introduced, since the radius of the ap- 
parent inner edge, Rf™ d derived by spectral fitting based 
on equation ([ll]) is generally overestimated by a factor of 
few because of negligence of the boundary term, so as to 
remove such artificial effects and to evaluate the radius of 
the real inner edge (Rf^ s = r]R™ od , see the fourth column 
of tables). In this paper, we set (£,7?) = (1.7,0.41) accord- 
ing to Kubota et al. (1998). The "7]° bs — Rf£ s " diagram 
are shown in Fig ^|. 

In the cases with \± = 0.01, the fitting results are nearly 
consistent those with p, = 0. As fi increases, i?? bs decreases 
(see Tables 1-3). For \i > 0.4, in particular, the inner 
disk temperature can never be larger than 2.0 keV. As 
we described in the previous section, when the effective 
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Fig. 4. The fitting parameters T ; ° ,R°^ B are plotted in this 
figure. The black-hole mass is M = IOMq . The solid lines are 
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Fig. 5. The thermal equilibrium (S-shaped) curves for dif- 
ferent parameter sets; (a, /i)=(1.0, 0.0), (0.1, 0.0), (0.1, 0.1), 
(0.1, 0.4), and (0.01, 0.0), respectively. The abscissa is the 
normalized mass accretion rate, and the ordinate is the sur- 
face density both on the logarithmic scale. The dashed-dotted 
lines represent the loci of H/R = 1 for several value of a. 

a is large, the disk cannot release its gravitational energy 
efficiently even for rh = 100. 

It is difficult to discriminate different models from the 
observational fitting results in Fig [|, because the differ- 
ence between fj, = and 0.4 is only a factor of a few. In 
the case of smaller rh value, there is practically no change. 
When rh is large, rh ^ 30, on the other hand, 7]° bs changes 
by factor of a few, and Rf^ s changes by one order of mag- 
nitude, depending on fi values. 

4. Discussion 

We thus expect different time-dependent behavior of a 
single source, depending on the ^values, if large rh modu- 
lation is observed. In the case of //~ 0.5, importantly, the 
disk inner radius Rf^ s does not change much even when rh 
changes a lot, in contradiction with the observed behav- 
ior of IC 342 which clearly shows a decrease in Rf^ s when 
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L increases (Watarai et al. 2001). Thus, the case with 
fj, 0.5 can be ruled out at least for the case of ULXs. 

According to the SSD theory, the a parameter affects 
the flow velocity (v r ) and thus matter density (£), in such 
a way that mass-flow rate (oc v r Ti) is kept constant, but 
does not largely affect its local emitting flux (as long as 
a simple blackbody radiation is assumed). In the high 
accretion-rate systems, however, the temperature within 
the transonic region of the flow is directly affected by the a 
value. Unlike the SSD the thickness of the high-luminosity 
disk becomes mildly thick (the relative thickness H/R is 
less than but of order unity) and Thomson optical depth 
is quite large. Therefore, photons generated by viscosity 
take longer time to reach the disk surface than the accre- 
tion timescale so that the generated energy by viscosity 
cannot be radiated immediately. 

From the observational point of view, the parame- 
ter acts more effectively on the fitting parameters in the 
slim disk than in the standard disk. This result leads the 
conclusion that the super-Eddington sources are useful for 
investigating the nature of viscosity. We find that, even 
if objects have truly high accretion rates, they may look 
like a standard-type disk, if fx is large. If we accurately 
measure the change in s with changes of m, we will be 
able to give a good estimate to the /i value. 

The origin of the viscosity in accretion disk is still poorly 
known, so it is practically impossible at the present to 
constrain the value of /z from the first principle. Also, 
the observational determination of /i is not easy. If u 
has a larger value in reality, the disk does not enter the 
advection-dominated regimes, even if the super-critical ac- 
cretion is realized. The results of the present study are 
useful for constraining the viscosity model. For example, 
large fi (> 0.1) models do not agree with the observed 
properties of ULXs, which show a tendency of decreasing 
R°n s (~ 3r g) at lar g e luminosity (Watarai et al. 2001). 
Hence, the viscosity model with a large /j parameter can 
be ruled out for such sources. 

Figure |^ illustrates the thermal equilibrium curves of 
the calculated disk at a fixed radius (r = 7r g ). The straight 
lines in the upper-left corner of the figure represent the 
loci of H/r = 1 for several values of a =0.01, 0.1, and 1.0 
from the bottom, respectively. How close the equilibrium 
curve is to the line of H/r = 1 is a good indicator to as- 
sess how the advection is important in energy equation, 
since Qadv ^ (H/r) 2 Q vis (Kato et al. 1998). The relative 
distance between the upper branch of the S-shaped curve 
and the line of H/r = 1 hardly changes for each set of the 
same a in Fig |5|. On the other hand, it is clearly demon- 
strated that when n is large, the kinks are not strongly 
manifested. Limit-cycle oscillations can occur at high m, 
only when \x is small, and the amplitudes of oscillations 
decreases as /i increases (Honma et al. 1991; Szuszkiewicz, 
Miller 1997, 1998). Large a values tend to stabilize the 
disk structure so that for even higher u(> 0.5) the limit 
cycle oscillations cannot occur. 

The occurrence of limit-cycle oscillations in real systems 
is still in open question. There is a hint of such oscillation 
in GRS1915+105 (Yamaoka et al. 2001), but we need 



careful data analysis in future work. 

The complementary work to the present study is made 
by Manmoto, Mineshige (2001) who consider fully rcl- 
ativistic treatments of the slim-disk problem but only 
for the canonical viscosity prescription (t rtp — — aptot) 
with a = 0.1. Main results are: (i) for the disk around 
a Schwarzschild black hole Rf^ s decreases from 3r g to 
~ 1.5r g as mass-accretion rate increases beyond a critical 
value M cr it ~ 30Le/c 2 ; (ii) for the disk around a Kerr hole 
(KBH), conversely, Rf^ s increases from < r g to ~ 1.5r g 
as mass-accretion rate increases because of significant self- 
shielding of the innermost hot region; (iii) at accretion rate 
above the critical value, disks around Schwarzschild and 
Kerr holes look quite similar; (iv) in both cases, the tem- 
perature distribution is somewhat flatter, txr -1 / 2 , at large 
luminosity. These results are basically in good agreement 
with the present study based on the pseudo-Newtonian 
treatment (except for the case of Kerr holes which are not 
considered here), but there exist some quantitative differ- 
ences. In the present study, hence, we are more concerned 
with the qualitative trend and thus it is unlikely that our 
main conclusions significantly change, even if fully rela- 
tivistic treatments expect that very small Rf^ s < r g may 
be an artifact of our pseudo-Newtonian treatments. 

For more realistic study along this line we also need 
fully 2D studies, taking into account the radiative transfer 
towards the vertical direction and various MHD processes. 
This is left as future work. Time dependent modeling of 
bursting behavior observation GRS1915+105 is another 
interesting subject in future work. 

Recently, King et al. (2001) proposed that large L 
and high TJ° bs observed in ULXs could be understand, 
if X-ray is moderately beamed. We suggest that spectral 
fitting with the p — free model will be able to test this 
hypothesis, since it only the thermal emission from the 
innermost region is beamed, the resultant spectrum will 
be of a singlc-tcmpcrature blackbody; i.e., large p will be 
obtained. Otherwise, we expect 0.5 < p < 0.75. We can 
thus judge it beaming is substantial in future observations 
(with, e.g., XMM Newton). 

5. Conclusions 

We have discussed the properties of the slim disk for 
different viscosity prescriptions and their soft X-ray sig- 
natures. Our conclusions are summarized as below. 

1. Some a dependence is found in the radial temper- 
ature distribution of the innermost region of high 
accretion-rate systems (m > 30); an increase in a 
tends to increase the temperature of the transonic 
flow, making a spectrum somewhat harder, although 
the change is not large and thus not easy to detect 
observationally. 

2. When the parameter [i takes a larger value (e.g. 
fj, = 0.5), the slim-disk nature (small Rf^ s < 3r g and 
flatter temperature profile, T e g oc r -1 / 2 ) is not ap- 
preciable even at high luminosity comparable to the 
Eddington luminosity. 
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3. The origin of the viscosity is still unknown, but ob- 
servations of high luminosity objects will be able to 
discriminate different viscosity prescriptions for ac- 
cretion disks. For example, large /x values are not 
incompatible with the behavior of some ULXs. The 
viscosity should thus be more dependent on total 
pressure, rather than gas pressure solely. 

We are grateful to the referee, E. Szuszkiewicz, for 
pointing out the initial of errors in the first version and for 
useful comments and discussion, which helped our mak- 
ing the revised version. We would also like to thank 
M. Takcuchi and T. Kawaguchi for discussion. This 
work was supported in part by the Grants-in Aid of the 
Ministry of Education, Science, Sports, and Culture of 
Japan (13640238, SM). 
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K. Watarai and S. Mineshige 

Table 1. Results of fitting with p-free model (0.2-10keV). (M = 10M e , fj, = 0.1, a 



M/(L E /c 2 ) 


T^obs 
in 


(keV) 


R™ od (km) 




P 


1 


0.56 (0.2-8 


.0 keV) 


245.4 


3.35 


0.75 


3 




0.75 


237.1 


3.24 


0.74 


10 




1.02 


211.0 


2.88 


0.72 


33 




1.60 


142.1 


1.78 


0.68 


100 




2.86 


56.7 


0.77 


0.61 


333 




3.37 


40.1 


0.55 


0.58 



Table 2. Results of fitting with p-free model (M = 10M Q , /z = 0.4, a = 0.1) 



M/(L E /c 2 ) 


T° bs (keV) 


R™ d (km) 




P 


1 


0.54 (0.2-5.0keV) 


249.6 


3.41 


0.75 


3 


0.75 


247.0 


3.37 


0.74 


10 


0.99 


233.5 


3.19 


0.73 


33 


1.39 


200.1 


2.73 


0.71 


100 


1.90 


138.1 


1.89 


0.65 


333 


2.92 


63.8 


0.87 


0.59 



Table 3. Results of fitting with p-free model (M = 10M Q , fi = 0.49, a = 0.1) 



M/(L E /c 2 ) 


T obs ( ke y) 


R™ d (km) 


nobs / 
^Mn / ' S 


P 


1 


0.54 (0.2-8.0keV) 


241.77 


3.36 


0.74 


3 


0.73 


266.24 


3.69 


0.76 


10 


0.90 


290.92 


4.04 


0.75 


33 


1.16 


296.76 


4.10 


0.73 


100 


1.38 


276.28 


3.83 


0.68 


333 


1.48 


266.66 


3.07 


0.62 



Table 4. Results of fitting with p-free model (M = 1OM , fj, = 0, a = 0.01) 



M/(L E /c 2 ) 


T^obs 
in 


(keV) 


i?f n od (km) 


< s Ag 


P 


1 


0.56 (o.2-t 


i.OkeV) 


245.4 


3.41 


0.75 


3 




0.73 


236.4 


3.28 


0.74 


10 




1.04 


207.0 


2.87 


0.71 


33 




1.51 


153.5 


2.13 


0.68 


100 




2.50 


68.64 


0.95 


0.61 


333 




2.94 


51.59 


0.72 


0.57 



Table 5. Results of fitting with p-free model (M = 10M Q , /i = 0, a= 1.0) 



M/(L E /c 2 ) 


T obs ( ke y) 


R™ d (km) 


nobs / 
"in /'g 


P 


1 


0.56 (0.2-8.0keV) 


244.2 


3.39 


0.75 


3 


0.75 


234.9 


3.26 


0.74 


10 


1.04 


204.0 


2.83 


0.71 


33 


3.09 


46.04 


0.64 


0.65 


100 


3.93 


36.82 


0.51 


0.62 


333 


4.18 


30.67 


0.43 


0.58 



